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Abstract 

Currently, the structure of the A(3872) meson is unknown. Different competing models 
of the cc exotic state A(3872) exist, including the possibilities that this state is either 
a mesonic molecule with dominating D^D*^ + c.c. composition, a ccqq tetraquark, or a 
cc-gluon hybrid state. It is expected that the A(3872) state is rather strongly coupled to 
the pp channel and, therefore, can be produced in pp and pA collisions at PANDA. We 
propose to test the hypothetical molecular structure of A(3872) by studying the D or D* 
stripping reactions on a nuclear residue. 
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1. Introduction 

The discovery of exotic cc mesons at B-factories and at the Tevatron stimulated interest 
to explore the possible existence of tetraquark and molecular meson states. The famous 
X(3872) state has been originally found by BELLE [l| as a peak in invariant 

mass spectrum from exclusive -A- K^7i~^7r~ J/p decays. Nowadays the existence of the 
A(3872) state and its quantum numbers = I"'"’'' are well established j^. In particular, 
radiative decays A(3872) -A J/pf, A(3872) -A 'ip'{2S)'y point to the positive C-parity 
of the A(3872). Probably the most intriguing feature is that the mass of the A(3872) 
is within 1 MeV the sum of the and D*^ meson masses. This prompted the popular 
conception of the A(3872) being a DD* + DD* molecule. 

To probe the molecular nature of the A(3872) structure has been difficult. So far, 
most theoretical calculations have been focused on the description of radiative and isospin- 
violating decays of the A(3872). For example, the A(3872) -A J/pf decay can be well 
understood within the DD* + c.c. molecular hypothesis [^. On the other hand, the mea¬ 
sured large branching fraction i?(A(3872) —)■ p'{2S)p/ -A J/pp = 3.4 ± 1.4 
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seems to disfavour the molecular structure and requires a significant pure cc admixture 
in the X{3872) j^. The theoretical predictions for the decay rates are, however, quite sen¬ 
sitive to the model details even within various approaches like charmonium or DD* + c.c. 
molecular models. 

In this letter we suggest to test the charm meson molecular hypothesis of the X(3872) 
structure in pA collisions at PANDA. Assuming that the X(3872) is coupled to the pp 
channel, we consider the stripping reaction of the D-meson on a nuclear target nucleon 
such that a D* is produced and vice versa. We show that the distribution of the pro¬ 
duced charmed meson in the light cone momentum fraction a with z-axis along p beam 
momentum, 

2{ujd*0^d*) + kf),) 

a = — - - , (1) 

Ep + mN+ Plab 

will be sharply peaked at a ~ 1 at small transverse momenta which allows to unam¬ 
biguously identify the weakly coupled DD* + c.c. molecule. Here, k^* and a;D*(k£)») = 
(k|,, -|- are, respectively, the momentum and energy of the produced D* meson in 

the target nucleus rest frame. Similar studies of hadron-, lepton-, and nucleus-deuteron 
interactions at high energy have been proposed long ago to test the deuteron structure 
at short distances as in the spectator kinematics the n- or p-stripping cross sections are 
proportional to the square of the deuteron wave function. For the X(3872) this idea is 
depicted in Fig. [U (details follow below). 

2. X(3872)-proton cross section 

For brevity, the bar, which can be seen over the D* or D, will be dropped in many 
cases below. The charge conjugated states are implicitly included in the calculated cross 
sections. 

The most important ingredients of our calculations are the total Xp cross section and 
the momentum differential cross section Xp —)■ D*(D)-|-anything. In the molecular picture, 
the latter cross section is the D(D*)-meson stripping cross section. To calculate the total 
Xp cross section within the Glauber theory, we start from the graphs shown in Fig. [1] 
which assume the DD* composition of X(3872). It is convenient to perform calculations 
in the DD* molecule center-of-mass (c.m.) frame with proton momentum pp directed 
along z-axis. The invariant forward scattering amplitudes of the hrst two processes are 



[ (fk |?7(k)|^Mp£,(0) , 

J ujd 

( 2 ) 


j d?k |V’(k)p*Mpo*(0) , 

J OJd* 

( 3 ) 


where mx = ooc + ujd* is the mass of the molecule and ujd {ujd*) is the energy of D 
(D*)-meson. (The different assumptions on the momentum dependence of meson energies 
discussed in the next section have practically no effect on the Xp cross section.) The 
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Figure 1: Processes contributing to the forward scattering amplitude of a proton on the DD* molecule. 
Wavy lines denote the pD and pD* elastic scattering amplitudes. Straight lines are labeled with particle’s 
four-momenta. The blobs represent the wave function of the molecule. 


molecule wave function in momentum space is defined as 

/ cl‘^T 


( 4 ) 


where k is the D* momentum in the DD* c.m. frame, with the normalization condition 

/ (i^/c|'0(k)p = 1. 

For the calculation of the third and forth processes in Fig. [T]we apply the generalized 
eikonal approximation (GEA) which assumes the nonrelativistic motion of D and 

D* inside the molecule. In this approximation, the propagator of the intermediate proton 
depends only on the .^-component of momentum transfer q = k^i* — k'^,*, while the pD and 
pD* elastic scattering amplitudes depend only on the momenta of incoming particles and 
on the transverse momentum transfer. Thus, we obtain 




(Pk(Pq imx 


( 27 r )^ 2ujdojd* 
Pkdp'q imx 
(27r)^ 2udOJd* 


2pp{q^ + ie) 


,MW(0)= /- W, 

o,. ,. 2pp{-q^ + ie) 


( 5 ) 

( 6 ) 


Therefore, 


iM^^\0) + iM^P (0) 


cPkcPqt 


- - 'ijj*(k - qt)iMpD*{qt)iMpD{-qt)'ip(k) . 

AujdUJd*Pp 


( 7 ) 
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The optical theorem for the proton-molecule forward scattering amplitude is 


ImM(O) = 2p^mx(T^°x ■ (8) 

Substituting M(0) = and using the parameterization 

of the strong interaction scattering amplitudes in the usual form as 

with Jp^(*) (k^ic*)) = [{EpCj^i*) — Ppfc^{*))^ — being the Moeller flux factor we 

obtain the following expression for the proton-molecule total cross section: 


<5 = j d^mk)ni,o(-k)a';i+i^o-(kKi.] -\j d:>ki,{k)i^o{-ky^i,o-{k)<j'^i. 


X 




( 10 ) 


where the normalized flux factors are dehned as Ip^(*)(k) = Jp^(*)(k)/ppCi;£)(,). In the 
small binding energy limit the molecule wave function decreases rapidly with increasing 
momentum k and becomes negligibly small at fc ^po'^- I'his case one can set 5^0 = 
BpD* = 0 and perform the Taylor expansion of the flux factors in in Eq. fllOp . Then, for 
the S-state molecule with accuracy up to the linear terms in and assuming that 

mo — Tn*j^, ~ we obtain the formula 


cr, 


tot 

pX 


= cr, 


tot 

pD* 


I ^tot 
+ ^pD - 


—tot ^tot 
^pD*^pD / -2\ 

\r )dd* 


4:71 


( 11 ) 


in line with previous calculations of the proton-deuteron total cross section 

We choose the wave function of a DD* molecule as the asymptotic solution of the 
Schroedinger equation at large distances: 


Tp^r) 


e 

V 27r r 


( 12 ) 


where the range parameter k = ^/2jJEb depends on the reduced mass p = rnDniD* / {'mo + 
mo*) and on the binding energy Eb of the molecule. The corresponding momentum space 
wave function is 

V'(k) = • (13) 

Let us now discuss the input parameters of our model. Since there is no experimental 
information on Dp and D*p interactions, we rely on simple estimates in the high-energy 
limit. For small-size qq conhgurations the color dipole model predicts the scaling of the total 
meson-nucleon cross section with the average square of the transverse distance between 
quark and antiquark in the meson, which is proportional to the square of the Bohr radius 
tb = 3/4pas. Here, p = mqmq/{niq -|- mq) is the reduced mass with mq and mq being the 
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constituent quark and antiquark masses. The Bohr radii of pion, kaon, Zi)-meson and J/'ij) 
are ordered as tb-k > tbk > ^bd > 'f'Bj/ip- Hence, we expect that the total meson-nucleon 
cross sections follow the same order. At a beam momentum of 3.5 GeV/c (1/2 of the 
momentum of A (3872) formed in the pp ^ X process on the proton at rest) the total 
71^p and K^p cross sections are about 28 mb and 17 mb, respectively j^. The J/pjp cross 
section is expected to be much smaller, 3.5 — 6 mb (c.f. |9| and refs, therein.). We assume 
the total Dp cross section = 14 mb, i.e. slightly below the K^p total cross section. 
This choice is in reasonable agreement with effective held theory calculations (loj |. 

It is well known that at incident energies of a few GeV, the amplitude of meson (nucleon) 
- nucleon elastic scattering is (to a good approximation) proportional to the product of the 
electric form factors of the colliding hadrons (see e.g. ll| and refs, therein). Thus, in the 
exponential approximation for the f-dependence of the form factors, the slope parameters 
BpM of the transverse momentum dependence of the meson-proton cross section at small t 
should be proportional to {r"^)p+{r'^) m i where {r‘^)p and are the mean-squared charge 

radii of the proton and meson, respectively. Since {r‘^)M slope parameters 

should be also ordered as the Bohr radii. Empirical values at piab = 3.65 GeV/c are 

and BpK+ = 4.12 ± 0.12 GeV“^ as htted at 0.05 < —t < 0.44 
On the other hand, Bpj/^ = 3 GeV“^ at the comparable beam momenta 0 . 
We will assume the value Bpr, = 4 GeV since the Bohr radii of kaon and H-meson 
differ by ~ 30% only. For the pD* interaction we assume for simplicity = o'/d 
BpD* BpB)‘ 


Bp^+ = 6.75 ± 0.12 GeV-2 


GeV2 12 


Our educated guess on the D- and Zl*-meson-nucleon cross sections and slope parame¬ 
ters should of course be checked experimentally. The empirical information on can be 
obtained by measuring the A-dependence of the transparency ratio of H-meson production 
in pA reactions at beam momenta beyond the charmonium resonance peaks, where the 
background pp —?■ DD channel dominates. The slope parameter Bp^ can be addressed by 
measuring the transverse momentum spread of H-meson production in pA reactions. 

We will further assume that the A(3872) wave function contains 86% of D^D*^ + c.c. 
and 12% of the D^D*~ + c.c. component as predicted by the local hidden gauge approach 
. The binding enercy of D^D*^ is likely less than 1 MeV [l^ and can not be determined 
from existing data [2| accurately enough. We set = 0.5 MeV and = 8 

MeV in numerical calculations. This corresponds to the range parameters Kbod*o = 0.16 
fm“^ and kd+d*- = 0.64 fm“^. With these parameters the total pX cross section (fTOj) 
is a^x — 26 and 23 mb for D^D*^ and D~^D*~ components, respectively, at the molecule 
momentum of 7 GeV/c in the proton rest frame. 


3. D(D*) stripping cross section 


In high energy hadron-deuteron reactions, the main contribution to the fast backward 
nucleon production (in the deuteron rest frame or equivalently - fast forward in the deuteron 
projectile case) is given by the inelastic interaction of the hadron with second nucleon of 
the deuteron [l^ . For large nucleon momenta the spectrum is modihed as compared to the 


impulse approximation (lA) due to the Glauber screening and antiscreening corrections 15 
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Figure 2: The amplitude for the process X(3872) + p ^ D* + F where F = {Fi,... ,Fn} is an arbitrary 
final state in the pD interaction. See also caption to Fig. [TJ 


since the hadron may interact with both nucleons. In a similar way, in calculations of the 
cross section X+p —)■ +anything, we take into account the lA diagram (Fig. [2^) and the 
single-rescattering diagrams of the incoming proton (Fig. [2 ]d) and of the outgoing proton 
or of the most energetic forward going baryon emerging from the inelastic pD interaction 
(Fig. |2t). The expressions for the invariant matrix elements for the processes (a) and (b) 
in Fig. [2] are straightforward to obtain in the c.m. frame of the molecule state X: 




D* 

Ud 
1/2 


zm 


X 


2ppy/%J]^Lj]p* 




d^rzjj{r)Q{—z) 


(qi) 


(14) 

(15) 


where k = k^*. In the case of we applied the GEA by expressing the propagator of 
the intermediate proton in the eikonal form and using the coordinate representation with 
r = rp,* — rp. The explicit form of the amplitude can be written only for specihc 
outgoing states F. However, for the diffractive states including the leading proton, the 
expression for can be obtained from the expression for by replacing 0(— 2 ) —)■ 
Q{z), which reflects the change of the time order of the pD* and pD interactions. Thus, 
for the diffractive outgoing state F the expression for is given by Eq. flT5|) 

with replacement 0(—z) 1 (neglecting small differences in momenta of incoming and 

outgoing proton in elementary amplitudes). We assume that the same replacement can be 
done for any hnal state F. By summing over all states F we then obtain the momentum 
differential D* production (i.e. D-stripping) cross section in the molecule rest frame: 


da, 


d^k 


pX—^D* 


D* 


{2Ti)^2uD*4:Ppmx 


^ J 

spins and sorts of 

d^VF 


{px -b k'p, -Pp-px) 


+ +M(")|2(27r)^ 

d^PF 


{2xf2Ex, {2xf2Ex^ ’ 


(16) 


where px is the four momentum of the molecule {p\ = m\). With a help of the unitarity 
relation for the elementary amplitudes [l^ the sum over spin states and sorts of F and the 
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integration over phase space volume can be reduced to the products of the imaginary parts 
of elastic scattering amplitudes. This leads to the following expression for the momentum 
differential cross section in the molecule rest frame: 


d^O'pX^D* 

¥k 


K 


<^X,Z5(-k)|^(k)|2«: , 


1 


- ap%IpD*(k) 


(27r)2 '^*(k) 

‘ d‘^qtd'^q[ '0(k + qt)t^*(k + q') 
{27tY |^(k)P 

Xe“['®p-D*(q?+qf )+SpD(ql-qt)h/2 


(17) 


(18) 


The hrst term in the r.h.s. of Eq. ffTSj) is the pure lA contribution. The second and 
third terms are, respectively, the screening and antiscreening corrections (see Eqs. (8a) 


and (8b) in [l5|). The D* meson is assumed to be on its vacuum mass shell, a;£)*(k) = 

\/'ml 


|- k^, while the energy of the D meson is calculated from energy conservation, 
IjJD (-k) = mx — a;D*(k). (The condition loo > 0 constrains the maximum momentum 
of the emitted D*, k < 3.3 GeV/c. Above this value our model looses its applicability.) 


^tot 


tot 


In the case of the X-meson production one has to exchange IpD IpD*, 0‘p'b -H- 
and BpD -H- Bpo* in Eqs. (IT7)) . (lT8|) . In this case the on-shell condition is applied to the 
D-meson, while the D* energy is determined by energy conservation. 

It is convenient to express the differential invariant D* production cross section flT7)) in 
terms of the relative fraction a of the light cone momentum of the DD* molecule carried 
by the D*: 


d?a, 


U)£)*- 


pX^D* 


d?a, 


d?k 


= a- 


pX—^D* 


dad?kf 


= Gy“‘(a,k,) , 


(19) 


where a = 2 ( 0 ;/)* (k) — k^)/mx- Figures |3] and |4] show the differential cross section of 
D*^ and production from X(3872) collisions at 7 GeV/c with proton at rest as a 
function of a for several values of transverse momentum kf. At kt = 0, the cross section 
has a sharp maximum at a ~ 2mj:,*/mx — 1.04 and is almost unaffected by the screening 
and antiscreening corrections. With increasing the width of a-distribution increases 
while the screening and antiscreening corrections to the lA term become important. This 
is expected since the large-fc* component of the molecule wave function corresponds to 
small transverse separation between D and D*. The corrections become large for a ~ 1 
and large transverse momenta as can be directly seen from the structure of the integrands 
in Eq. ffTSj) . Indeed, a ~ 1 corresponds to /c^ ~ 0 in the molecule rest frame. Then at 
hnite transverse momentum transfer qt the ratio '^^(kj -|- /'ip* {k.t) is less than unity at 

/ct = 0 and asymptotically tends to unity with growing kt- Due to the extremely narrow 
wave function of the D^D*^ molecule in momentum space, the screening and antiscreening 
corrections are sharply peaked at a ~ 1.1 and develop structures in the a-dependence of 
the cross section at large transverse momenta. In the case of pA reactions these structures 
are slightly smeared out due to the nucleon Fermi motion (see Fig. [6] below). 
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X(7 GeV/c) + p ^ D*° 
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Figure 3: The invariant differential cross section of D*^ production in X(3872)p collisions at puh = 7 
GeV/c. Thick solid line - full calculation according to Eas. (|T7)) - (fT^ . Thin solid line - the calculation 
taking into account only IA and screening term of Eo. lfT^ . Dashed line - the calculation with k = 1 in 
Eq.dni, i.e. only with the lA term. The inset at /cj = 0 shows the behaviour of the differential cross 
section for a smaller range of a. 
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Figure 4: Same as Fig.[3]but for D*^ production. 
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Figure 5: The amplitude of the process pA —>■ D*T{A — 2)*. The wave functions of the initial and final 
nuclei are denoted as ijjA and 'ipA-2, respectively. The mass numbers are shown as subscripts. Wavy lines 
represent elastic scattering amplitudes on nucleons. stands for the arbitrary final state particles in 
the semi-inclusive process X2 —>■ D*J-. The summation is performed over all possible sets of nucleon 
scatterers {ni},{n 2 } and {na} for the p, X and D*, respectively. 


4. D* and D production off nucleus 

In antiproton-nucleus interactions, we focus on the D* (or D) meson production in 
the two-step process pp -A X, XN -A D*{D) + anything. Similar to the case of Xp 
interactions, we apply the Glauber theory to calculate the differential cross sections of the 
D* production in antiproton-nucleus interactions. We start from the multiple scattering 
diagram shown in Fig. |5] which can be evaluated within the GEA. We will assume that 
the nucleus can be described within the independent particle model disregarding the c.m. 
motion corrections (c.f. 0 ). The incoming antiproton, intermediate molecular state X 
and outgoing D*-meson are allowed to rescatter on nucleons elastically an arbitrary number 
of times. The D* production cross section is proportional to the product of the sum of the 
amplitudes of Fig. |5]and their conjugated. The X state is formed on a proton 1, while the 
D* is produced in the collision of X with a nucleon 2. The nucleons 1 and 2 are hxed in the 
direct and conjugated amplitudes while the sets of other nucleon scatterers are arbitrary. 
The leading order contribution is given by the product term without elastic rescatterings. 
Nuclear absorption corrections are accounted for by summing all possible product terms 
with non-overlapping sets of nucleon scatterers. This gives the following expression for the 
momentum differential cross section of D* production on the nucleus: 

/ p d^FW'^ * n 

dViPp,surv(bi, -oo, zi) J kt - -pu) 

OO 

X j dZ2Vx,smv{i>l,Zl,Z2)pihi,Z2)VD*,surv{hi,Z2,00) , (20) 

21 


^d^apA^D* 

dad^kt 
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where 


(Pplt 


\Mx-pi\^Vp 


np(ri;pii, 


( 21 ) 


is the in-medium width of p with respect to production of X with transverse momentum 
Pn; Vp = P\ah/Ep is the antiproton velocity; np(ri; pii, is the proton occupation 

number; Ei = mx — B with B = 8.6 MeV being the nucleon binding energy in ^°Ar 
nucleus. The longitudinal momentum of the proton 1 is obtained from the condition 
of on-shell production of the state X in the process p 1 —)■ X: 



+ El + 2EpEi - m\ 
2piab 


( 22 ) 


The nucleon occupation numbers are taken as the depleted Fermi distributions supple¬ 
mented^^ high-momentum tail due to short-range quasideuteron correlations (SRCs) 

300: 


nq(r;p) 


(1 - P2,q)Q{pF,q - p) + 


'^^P2,qPq\^d{p) P0(P - PF,q) 

OO ’ 

/ dp'p'‘^\pJd{p'W 


PF,q 


q = p,n (23) 


where pF^q{r) = [ 37 r^pq(r)]^A the nucleon Fermi momenta, P 2 ,p = 0.25 and P 2 ,n = 
P 2 ^pZ/N are the proton and neutron fractions above the Fermi surface, Pq{r) are the nucleon 
densities, and iJd{p) is the deuteron wave function. In Eq. fl20|) . the nuclear absorption is 
given by the survival probabilities of the antiproton, the molecule, and the D*\ 


'Rp,surv(bl, CX), 2^i) 


P XjSurv (bl , 2^1 , ^ 2 ) 


PD*,surv (bi ,Z2,0C) 


exp 


exp 


-a/p / dzp(hi,z) , 


— OO 

Z2 


-a/x / dzp(buz) )- , 


•21 


exp < / dzp{hi,z) > , 


•22 


(24) 

(25) 

(26) 


where p = Pp -|- is the total nucleon density. We use the two-parameter Fermi distribu¬ 
tions of protons and neutrons j^. As usual in the Glauber theory, Eqs. fl2T|) - fl26|) neglect 
the Fermi motion of nucleon scatterers. In a similar way, in writing Eq. fl20l) we neglected 
the Fermi motion of nucleon 2 since the elementary cross section flT^ depends only weakly 
on the proton momentum (via the flux factors, screening- and antiscreening contributions) 
and is in leading order proportional to the square of the molecule wave function. However, 
the transverse Fermi motion of proton 1 is taken into account in Eq. fl20l) in the high-energy 
approximation (c.f. 0). The latter implies that the light cone momentum fraction a can 
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be expressed in the target nucleus rest frame according to Eq. ([T]) where the Fermi motion of 
the proton 1 is still neglected. (We have numerically checked that using the exact Lorentz 
transformation to the c.m. frame of X to evaluate the invariant cross section ud* 
instead of using the inhnite momentum frame in Eq. fl20|) which conserves a and assumes 
Galilean transformation for kj produces indistinguishable results.) 

The pp ^ X matrix element in Eq. fl2ip is one of the major uncertainties in our calcu¬ 
lations. Its modulus squared can be formally expressed in terms of the partial decay width 
T'x—^pp 3-8 


\Mx;pl\^ = 


47r(2J 


X 


l)mlT 


X^ X^pp 




4m% 


(27) 


where the overline means averaging over antiproton and proton helicities and summation 
over the helicity of X. There is no experimental data on the partial decay width ^)^ 3872 )^pp- 
(The recent LHCb data on pp invariant mass spectra from B~^ ppK^ decays [211 do not 
allow to clearly identify X(3872) in the pp decay channel due to statistical limitations.) 
In the present calculations, we will use the value rx( 3872 )^pp — 30 eV as suggested by 
theoretical estimates j^. This value is about two times smaller than However, 

one should note that, in the molecular picture, the decay of the X(3872) to the pp state 
requires the production of only two qq pairs, and not three qq pairs as in the ordinary 
charmonium decay to the pp channel. Thus, the partial decay width of the X(3872) into 


the pp channel may be even larger than that of the Xci{^P) state 22 


Formula fl20|) has a simple physical interpretation if we express the integral j d^ri 
as J (Pbi J dzi- The factor Pp,gurv(bi, —cxo, zi) is the probability that the incoming from 
2 = —cx) antiproton with impact parameter bi will reach the point 2 ; = 2 ;i. The combina¬ 
tion [dzi/vp)d?pitd?V^^/d?pit is the W(3872) formation probability within the transverse 
momentum element d^pu when the p is passing the longitudinal element dzi. The factor 
’Px,surv(bi, 2 : 1 , 2 ^ 2 ) is the probability that the molecule will reach the point z = ^ 2 - The 
combination dz 2 {da/a)d^ktG^^^ (a,kt — f Pit)p(bi, 2 : 2 ) is the probability that a D* will 
be produced in the kinematical element dad^kt when the X(3872) is passing the longi¬ 
tudinal element dz 2 - Finally, the factor T’D*,surv(bi, 2 ^ 2 , cxo) is the probability that the D* 
will escape from the nucleus. In the spirit of the eikonal approach, all particles propagate 
parallel to the beam direction. (For example, we assumed that the transverse momentum 
of the molecule, pit, does not influence its trajectory.) The integration over Z 2 can be taken 
with the explicit forms of the survival probabilities Eqs. fl2^ . (l2^ . As a result Eq. fl20|) takes 
the following simple form: 




a- 


pA-^D* 


dad?kt - al%) 


d suj-v(bi, CX), 2:1) / d pit~ 




d^p 


It 


xG^“^^‘(a,ki - ^Pit)[PD * ,surv (bi,Zi,cx) T*x,surv(bi, . 2 i, cx)] . (28) 


In Fig. in] we display the differential cross sections of charmed meson production in 
antiproton collisions with argon nucleus at 7 GeV/c. The D and D* cross sections are 
peaked at a ~ 2mD/mx = 0.96 and a ~ 2m]:,*/mx = 1.04, respectively, and behave in 
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Figure 6: The invariant differential cross sections of D*^ and production in p'^^Ar collisions 

at piab = 7 GeV/c. The calculations are done using full cross section A(3872)p —> D*{D) in Ea. (l28l) 
including the lA term as well as screening and antiscreening corrections (see Eqs. (HZli-dlll)). For kt = 0, 
the cross sections of D*^ and production are divided by a factor of 100. 
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Figure 7: The ratio of D* production cross sections in and X{3872)p collisions at piab = 7 GeV/c 

for several values of transverse momentum kt as a function of light cone momentum fraction a. The ratio 
is normalized at unity for a = 0.5. Panel (a) - D*^. Panel (b) - D*^. 


similar way as a function of a and kt- The widths of a-dependence of the D*^ and 
cross sections are much smaller and the peak values are much larger as compared to the 
and cross sections. The a-dependence of D* and D production in pA collisions is 
dominated by the elementary cross section (c.f. Figs. [3|4|1 . However, a closer look reveals 
significant differences between D* production on a nucleus and on a proton due to the 
Fermi motion. These are better visible in the ratio of the two cross sections depicted in 
Fig. [3 At fci = 0 the ratio has a minimum at a ~ 2m£)*/mx because in this case the 
contribution from target protons with hnite transverse momentum pu is suppressed by 
the factor |'^(ki — ^pi4)p/|'0(kj)p. However, with increasing kt this factor becomes larger 
than unity for comoving proton 1. This leads to the observed local maximum in the a- 
dependence for ~ 0.1 — 0.4 GeV/c. At large kt or for large deviations of a from unity 
the ratio tends to the constant value. 


5. Uncertainty and backgronnd 

The uncertainty of our calculations can be read from Fig. |8l The unknown binding en¬ 
ergy of the D^D*^ molecule is the main source of uncertainty as a weaker binding produces 
a narrower a-distribution and vice versa. As a consequence of the partial cancellation be¬ 
tween the survival probabilities of D* and molecule the reduction of the pD and pD* cross 
sections from 14 mb to 7 mb leads to the reduction of the peak of D*^ production cross 
section by ~ 15% only. Of course, on the top of these two effects there is an uncertainty 
due to the experimentally unknown width Tx^pp which enters the cross section (|28|) as an 
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Figure 8: The a-dependence of production at fcj = 0 in p'^^Ar collisions at piab = 7 GeV/c. Panel 
(a) shows calculations with default cross section = 14 mb. The signal cross section (l28l) is shown for 
the different binding energies Ej, of the molecule. The background cross section (1291) is shown with 

SRCs (default calculation) and without SRCs. Panel (b) shows calculations with default Et = 0.5 MeV 
for the different ap°^ as indicated. It is always assumed that . Insets show the narrower region 

of a. The background cross section is divided by a factor of 3. 
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overall multiplication factor. The qualitative behavior of the a-distribution is not changed 
by varying the model parameters. 

The major background is given by the direct process pN DD*,DD* on the bound 
nucleon, because the thresholds of X(3872) and DD* production in pp collisions are almost 
the same. The background cross section can be calculated as (derivation is similar to that 
ofEq. ([28])) 


^ f j2 dapi^D*D{£) 

dad?kt (27r)3Eipiabfcf) i dVt 

X j dViPp,surv(bi,-oo,Zi)?7,i(ri;pit, A°) 

^ D* ,surv (bi,Zi,CX))[l -PD,surv(bl,Zi,CX))] , (29) 


where 


£ = e(a, k( = 



/2[(2 — Q!)m|)* + am|, + 2k'^] 
\ a{2 — a) 


1/2 


(30) 


is the c.m. energy of D and D*; qpi{e) = (e^/4 — and qDD*{^) = [(^^ + rn\) — 

rn?jpt,Y/Ae^ ~ are the c.m. momenta of the colliding pN pair and of the produced 

DD* pair, respectively. The longitudinal momentum kp, of i7-meson in the nucleus rest 
frame is calculated from relation 


2{ujd + kp)) 

Ep + mN+ Plab ’ 


(31) 


with ud = y kpi^ + rajp- The longitudinal momentum of the nucleon, A^, is given by 
Eq. fl22l) with mx replaced by e. 

For the D*^ production on the proton, the near-threshold S'-wave cross section is 
app^o*D{£) = 2app^jj*ojjo(e), where the direct (non-resonant) cross section Cpp^jp^opjo^e) 
has been taken from ref. |22|] which is the only estimate of the discussed cross section avail¬ 
able in the literature. (We included the factor of 2 since our D*^ includes both physical 
states, D*^ and D*^.) The estimate of was obtained by using dimensional counting 
considerations to express the cross section of pp —?■ D*^D^ at high energies in terms of 
the cross section of pp —?■ K*~K^ which is known in the limited energy range. As the 
next step, the pp —?■ D*^D^ cross section was extrapolated in towards the threshold 
and multiplied by the S'-wave fraction /i=o — 9.3% which is regarded by the author of 
|22|himself as “a crude extrapolation”. Thus, we feel that the near-threshold estimate 
of j22| can be considered as an order of magnitude estimate. In the case of the D*^ pro¬ 
duction on the neutron, the only possible channel is pn —)■ D*^D~. Thus, we assume 

0'pn^D*D{.£) = (^Pp^D*0D°{^)- 

The result of calculation using Eq. fl29p is shown in Fig. [HI The dependence of the 
background cross section on the pD* and pD cross sections is quite modest and follows 
the tendency of the signal cross section. Thus, at fcj = 0, the sharp peak of D* production 
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at « = 1.04 due to the stripping reaction is clearly visible on the smooth background. 
The peak is almost not influenced by intramolecular screening and antiscreening effects. 
Moreover, we expect that the elastic rescattering of antiproton and produced particles on 
the nucleons will practically not change the D* and D spectra at small kt 0. Finally, the 
influence of the SRCs - which are always included in default calculations - is mainly in 
the reduction of the nucleon occupancies at small momenta. This results in 15% reduction 
of the background at a ~ 1 (Fig. [8^) and the correspoding rescaling of the signal cross 
section (not shown). 

The signal-to-background ratio at the peak stays almost constant as the mass number 
of the target nucleus varies between 20 and 208. The shape of the a-dependence of the 
signal cross section practically does not vary with mass number in that region. The mass 
dependence of the total D*^ production cross section due to the stripping reaction can be 
well approximated by formula ao*o = 16 pb • With the high luminosity mode at 

PANDA, L = 2- 10^^ cm“^ s“^, the estimated production rate due to the stripping reaction 
is about 60 D*^ events per hour. 


6. Discussion and conclusions 


We proposed the idea of D{D*) stripping from the X(3872) state, to investigate if 
W(3872) has a molecular structure from the narrow peaks in a-distribution of D* and 
D at a ~ 1. Other microscopic models of W(3872), e.g. tetraquark or cc-gluon hybrid, 
would lead to the flat a-spectrum of D*{D). In such models, there are no primordial 
hadronic components in W(3872) to be “released” or “knocked-out”. Thus, the momentum 
distribution of D*{D) in the process XN ^ D*{D) would be dominated by phase space 
of the final state particles. 

There is, however, another possible source of narrow peaks in a-distributions of D* and 
D. The BELLE collaboration has found a significant near-threshold enhancement in 
the D*^I)^ invariant mass spectrum from B D*^D^K decays. We note that this does not 
exclude the existence of the D*^D^ bound state. (One similar example is A(1405) which 
lies about 30 MeV below K~p threshold and can be treated as a K~p quasibound state 
although it strongly influences the K~p —)■ and K~p —)■ cross sections at small 


beam momenta 2J, l25|].) But it is also possible that X(3872) is a resonance coupled to 


the D*^D^ + c.c. channel. If such a resonance state is produced in peripheral pA collisions, 
it will decay far away from the nucleus, since the width of X(3872) is less than 1 MeV. 
The resulting a-distributions of D*^ and will be also sharply peaked near a ~ 1 at 
small kt- However, in this case both decay products can in principle be detected. This 
gives a clear experimental signature for distinguishing such decay events. In contrast, the 
stripping events would contain only one meson, D*^ or Z)°, in the same kinematical region. 

The stripping reaction can be also considered for other production channels of X(3872), 
e.g. in proton-, electron- and photon-induced reactions on nuclei. Other exotic X,Y,Z 
states, such as the X(3940) [^, V(4140) j^, X(4160) (c.f. recent reviews j^, 13 

for a more complete list), may be interpreted as molecular states of D*D* or D*D*. These 
hypothetical molecular structures may also be tested by using stripping reactions, similar to 
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X(3872). Experimentally, such studies could be performed at J-PARC, FAIR, SPS@CERN, 
and EIC. 

Apart from cc exotic states, there are other mesons which possibly have molecular 
structures. The ao(980) and /o(980) states viewed as KK bound state can be produced in 
-y{7i)N fN (/ = oo, /o) reaction on the bound nucleon followed by the stripping process 
fN —)■ K{K) + anything on another nucleon of the nuclear target residue. The D*q{2317) 
viewed as DK and i7si(2460) viewed as D*K can be produced in pp DfD*^ (^?i) 
reaction on the bound proton followed by the stripping process —)■ D{D*) + 

anything or K + anything on another nucleon. The antiproton-nucleus 

interactions open another unique opportunity to produce fast antibaryons. In this way, 
e.g. the molecular K^p hypothesis for the A(1405) can be tested by using the two-step 
processes pp A(1405)A, A(1405)A^ —)■ -\- anything. Such kind of processes can be 

studied at PANDA and possibly at J-PARC. 

In conclusion, we have demonstrated that the spectra of D* and D in the light cone 
momentum fraction at small transverse momenta allow to test the hypothetical DD* molec¬ 
ular structure of the X(3872) produced in pA collisions at threshold. We propose to search 
the narrow peak in D* or D production at a ~ 1 and small kt as an unambiguous signal 
of the DD* molecular state formation in pA collisions in PANDA experiment at FAIR. 
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